Lorentz-Lorenz shift in a Bose-Einstein condensate 
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We study the quantum field theory of light-matter interactions for quantum degenerate atomic 
gases at low light intensity. We argue that the contact interactions between atoms emerging in the 
dipole gauge may be ignored. Specifically, they are canceled by concurrent infinite level shifts of the 
atoms. Our development yields the classic Lorentz-Lorenz local-field shift of the atomic resonance. 
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' The understanding of the interactions of light with matter in quantum degenerate systems has become especially 
ON \ topical after the first evidence for a Bose-Einstein condensate (BEC) of an atomic gas has appeared [Q. Recently, 
~~ * ■ Andrews et al. g| have also reported on non-destructive optical detection of a Bose condensate. 

In the limit of large detuning of the driving light from atomic resonance, the dynamics of the light and matter 
fields may be decoupled. It then turns out that the spectrum of the scattered light conveys direct signatures of atom 
statistics pH] , and that under various conditions even the phase of the macroscopic wave function may be observed 
optically [5|. On the other hand, for a dense enough sample and small enough atom- field detuning, the analysis of the 
CO ' response of matter requires a concurrent treatment of light with its own dynamics. The result is that nearby atoms 
alter the optical response of each other. A large body of work in this direction j^jlj] seems to have brought about 
thcgeneral notion that atom statistics should have only a minor effect on absorption, dispersion or diffraction of light 
[@,[l3j. Nonetheless, even for the Maxwell-Boltzmann gas and for the simplest optical properties such as refractive 
index, there still are no proven solutions for microscopic theories regarding a dense, near-resonance sample. 

In this paper we continue our rigorous quantum field theoretical analysis of light-matter interactions PJT^]. We 
_ , argue that the contact interaction terms between different atoms that arise in the length gauge do not have any effect 
on light-matter dynamics. As a result, for a model BEC we find precisely the Lorentz-Lorenz (LL) shift familiar from 
r*"** ! classical electrodynamics |l4]-|T6| . 

In the many-particle formalism of light-matter interactions, the Hamiltonian is often transformed into the length 
gauge using the Power-Zienau-Woolley transformation j[7). The electric displacement is then the basic dynamical 
degree of freedom, instead of the electric field, and the Hamiltonian picks up a polarization self-energy term 
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where P(r) = Yli^i^\ r i ~ r ) ^ s the electric polarization in the dipole approximation and d; and r$ are the dipole 
operator and the center-of-mass (cm.) position operator for the i th atom. On the right hand side of Eq. ([!]) we 
have ignored the divergent self-energies with i = j. Polarization self-energy describes a contact interaction between 
polarized atoms. One might argue that for a quantum degenerate Bose gas this term is important, as the de Broglie 
wavelengths of the atoms are of the order of the interatomic distances and the wave functions of different atoms overlap 
strongly. The polarization energy has been duly included in many theoretical approaches to light-atom interactions 
in the quantum degenerate regime [p| JTl|JT^ , |l8t . 

In our previous paper |l2] | we have derived a hierarchy of equations of motion for correlation functions that contain 
one excited-atom field and one, three, five, etc., ground state atom fields, for the limit of low light intensity. All 
contact interactions, such as those produced by the polarization self-energy (|l|), were included. For the total electric 
field E + we found the following monochromatic expression in the presence of scattering by atomic dipole moments 

e E+(r) = D+(r) + - / dV G(r - r') P+(r') , (2) 
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where is the positive frequency component of the driving electric displacement, the scalar constant k is defined in 
terms of the reduced dipole moment D as k = T> 2 /heo, and the positive frequency part of the polarization operator 
P + is given in second quantization by 

P+(r)=d 9e V4(r)^ e (r). (3) 

Here vb g and ip e are the ground state and the excited state field operators in the Heisenberg picture, and d ge is the 
dipole matrix element for the transition g — » e. In our notation the same level index appearing twice in a product 
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indicates a summation over the magnetic substates of the level. The monochromatic version of the 3x3 tensor 
propagator G(r) is 



Gij(r) = in 



dri drj 



e ikr 
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The expression G(r — r')T> is equal to the positive- frequency component of the electric field from a monochromatic 
dipole with the complex amplitude X>, given that the dipole resides at r' and the field is observed at r |jl6). The 
explicit expression is 

G(r)2> = ^{i 2 (nxD)xn^ + [3n(n • V) - T>\ (1 - %)^} -^VS(r), (5) 

with n = r/r and k — f2/c, Q being the dominant frequency of the incident light field. The volume integral over 1/r 3 
in Eq. (^) is not absolutely convergent in the neighborhood of the origin. The expression (|^) should be understood in 
such a way that the integral of the term inside the braces over an infinitesimal volume enclosing the origin vanishes [ p"2"| . 
Equation (Q) is precisely the classical expression of the dipolar field, including the peculiar delta function divergence 
at the origin |Tq |. 
With the definitions 

P,(n, • • ■ , r,_i; r,) = (^( ri ) . . . ^(r^OP+fa^to-i) . . . (6a) 
ptin, ■ ■ ■ ,r,) = (^( ri ) . . . v>t(r,)^( ri ) • ■ • ^(n)>, (6b) 

and considering a J g = — ► J e = 1 transition for simplicity, we were able to cast the hierarchy of equations for 
atomic correlations functions into the following form |l2| ] 

Pi(ri,---,r ; _i;r z ) = (iS - 7)Pj(r l5 ■ • ■ , rj_ x ; rj) 
l-l 

+ ^2 G(r t - r fe )P;(ri, • • • , r k -i, r fc+ i, • • • , r z ; r fc ) + inpi(r lr ■ ■ , r ; )Dj(r ; ) 

fc=i 

+ J d 3 ri +1 Q(ri - r ;+1 )P ;+1 (r 1; • ■ • , r ; ; r i+1 ) . (7) 

The quantity P; reflects correlations between the dipole moment of one atom and the positions of I — 1 other atoms, 
and pi is simply the density correlation function for I ground state atoms. 

In Eq. (0) the atoms are regarded as point dipoles. As we have already emphasized in Ref. JT2| , this assumption 
must fail at short distances. Real atomic interaction potentials are thought to have a hard core, which prevents the 
atoms from overlapping. Mathematically, one expects that all atomic correlation functions such as P; and pi vanish 
if any two position arguments are the same. Evidently, the contact interaction should be omitted as inconsequential 
in Eq. (Q) . In view of Eq. (0) , this is done by removing the delta function contributions from the field propagator G 
by the substitution Gy — > G^- , with 

Gj J .(r) = G y (r)+t«<y^(r)/3. (8) 

The key mathematical insight of this paper is that the same substitution applies even for point dipoles. In other 
words, the resulting light-matter dynamics is exactly the same, whether we use the field propagator G or G' in Eq. (Q). 
To demonstrate this contention, we consider the steady-state solutions of Eqs. (0). For each I we have a system of 
I coupled equations for P/(- • ■ ; ri), P;(- • ■ ; r^),..., P;(- ■ • ; r;). The coupling arises from the collective line shifts and 
linewidths, those G terms in Eq. (0) that appear outside the integral. The idea is that, if we solve this system of 
equations for a fixed I, in P;(- • • ; r;) and indeed even in the seemingly more divergent expression G(rfc — r;) P;(- • ■ ; r;) 
with k ^ I the delta functions of the form S(jj — r s ) (j ^ s) cancel. The S functions are inconsequential in integrals 
of the type J d 3 n_i G(r;_i — r;) P;(- • • ; r;), and hence in the solution of the hierarchy ([?]) for Pi. 

We begin the detailed discussion by noting that in the derivation of Eqs. (j^) and (Q) in Ref. |Q we used a cutoff in 
the wave numbers q of the photons. Here we again invoke the high-frequency cutoff. The S function is then regarded 
an ordinary function, albeit one that is sharply peaked over a small distance scale a. Seemingly unorthodox operations 
such as divisions by a S function are therefore well defined. After the calculations we set a — > 0. 

The steady-state solutions for P;(- • • ; r^), k = 1, • • • , I, of Eq. (Q) are obtained from a system of linear equations 
Ax = b, where the I x I matrix A and the / x 1 vectors x and b are defined by 
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A kq = Skq - (1 - S kq ) — G(r fe - r 9 ); x k = P ; (- ■ ■ ;r fe ); 
in 

h = ap;( r i: ' ' ■ : r O D F( r fe) + — d 3 ri +1 G(r fe - r i+1 )P ;+ i(ri, • ■ • ,r ; ;r ;+ i) . 

IK J 

Here a = — k/ (6 + i-f) is the polarizability of a single atom. The solution for this linear system may be expressed using 
Cramer's rule Q. In particular, xi = P/(- • • ; r;) = dct [A^]/ det (A), where = A kq , for q 7^ Z, and A^ = bi. For 
any two indices j ^ s, the matrix A contains two entries G(vj — r s ), one on the row j and column s and the other 
on row s, column j. The determinant det (A) therefore is a quadratic polynomial of the propagator G(rj — r s ). The 
coefficient of G(r^ — r s ) 2 is a function of the position variables {r^}, and may technically be zero for some sets of 
coordinates {r.;}. In our argument we assume that such accidentally zeros, if any, are inconsequential in the physics. 
We thus posit that det (A) contains a term proportional to G(jCj — r s ) 2 , and hence is a proper quadratic polynomial 
of 5(rj — r s ). By the same token, det [A^] is a second order polynomial of 5(rj — r s ) for j ^ s, j / I, s ^ I, and a 
first order polynomial in 6(rj — r;), j ^ I. 

Let us first consider the case I = N, where N is the total number of atoms. All the expectation values in 
Eq. (||) vanish for / > N. Thus, Pat+i = and bj = apN~Dp(rj), for j = 1, • • • , N. By the preceding argument, 
Pn — det[A( N ^]/ det(A) is a rational expression whose denominator is a quadratic polynomial in each S(rj— r s ), j =/= s, 
whereas the numerator is a quadratic expression of 5(jCj — r s ) with j ^ s, j 5^ N , s ^= N, and a first-order polynomial 
in S(rj — rjy), j ^ N. In the expression G(rAr — rjv-i)-Pjv(ri, ■ • ■ , rjv_i;rjv) all the delta function divergences in the 
numerator are therefore canceled by equal (or higher) powers of the same delta functions in the denominator. 

Next consider the case I = N — 1. Because the expression G(r^ — rAr_i)P7v(. . . does not contain uncanceled 
delta functions, its integral over needed to calculate the (N — 1) x 1 vector b for / = N — 1 does not contain 
any delta function divergences either. With similar arguments as before we then conclude that also G(ryv-2 — 
rjv_i) Pjy_x(- ■ • ; r^v-i) is free from uncanceled delta functions. By repeating the procedure for all I we confirm that 
the expressions G(r;_i — r;) P/(- ■ • ; rj), for any I, do not contain delta functions which could affect the integrations. 

Because our derivation is very technical, we give a simple example for the case of two atoms. Then, P3 = 0. The 
ground state density and the polarization at ri and r2 are coupled by the collective linewidths and line shifts leading 
to the resonant dipole-dipole interaction [Q. For / = 2, An = A22 = 1, ^12 = A 2 \ = — aG(ri — r 2 )/iK, and P 2 (ri; r 2 ) 
may be expressed using Cramer's rule. In particular, 



J d 3 r 2 G(ri - r 2 ) P 2 (ri; r 2 ) 



.3 n \ 1 sD+(r 2 ) + aG(ri -r 2 )D+(ri)/iK 

d J r 2 G (n -r 2 )ap 2 (n,r 2 — h — r ~ r r-^ . (9 

f - (aG(ri - Y 2 )/iK,y 

With r 2 7^ ri, the delta functions in G(ri — r 2 ) vanish, and at r 2 = ri we may cancel [S(rx —r 2 )] 2 . The delta functions 
do not affect the integral because the highest power of 6(r\ —r 2 ) is two, in both the denominator and the numerator. 
Thus, the contact interactions do not affect the polarization Pi(ri), which is obtained by inserting the expression (||) 
into Eq. (0). 

We have shown that all delta functions contained in G cancel in Eq. (f7j) . A corollary of the derivation is that the 
polarization self-energy term (Q) does not have any effect on the light-matter dynamics for point dipoles, and may be 
ignored. 

There is a concurring physical explanation for our result. Comparison of Eqs. (^|) and (Q) shows that the delta 
function lumps with the detuning <5, and gives a divergent frequency shift. As the dipole at r k draws nearer to the 
dipole at r g and the electric fields of the dipoles on each other grow stronger, at the same time the dipoles are also 
shifted further and further away from resonance with one another. Within our present viewpoint it is the level shift 
due to resonant dipole-dipole interactions that dynamically curbs the interaction energy between the dipoles. 

After arguing that the delta functions may be ignored in Eq. (R), we now rewrite the example presented in Ref. p^ ] 
for the optical response of a homogeneous condensate with the contact interactions omitted. The condensate is 
assumed to fill the half-space z > with the constant density p; for an ideal condensate the density p immediately 
implies the other correlation functions as pi = p l . The optical response is solved in an approximation that ignores 
the collective line shifts and linewidths. Specifically, we replace G in our earlier argument by G', and then ignore 
all tensors G' that reside outside integrals in Eq. (]?]). This approximation is only guaranteed to be valid for dilute 
condensates, pA 3 -C 1, where A is the wavelength of light pofl . What the exact limit of validity is and what types of 
corrections emerge at high densities is at the moment largely unknown. 

We consider the steady-state solution of Eqs. (Q) with I = I,.... The initial free field is written D^(r) = 
Dp eexp (ikz). The following set of damped plane waves solves Eqs. (^): 
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l_ 0, otherwise, 

with k' (S(fc') > 0) and P being the variable parameters. This may be seen by evaluating the following integral for 

21 > 

Ak' z<2 



f d 3 r 2 e*-G'(r 1 -r 2 )-ee J 

Jz 2 >0 

(11) 



2fc 2 + k' 2 ik > zi k' 2 ikzi 

3(k' 2 -k 2 ) 2k(k-k>) 



The integrals on the right-hand side of (0) produce sums of two exponentials, the vacuum component oc e lfcz , which 
should cancel the free-field terms, and a remaining e lk z term, which pairs up with the correlation functions P;. All 
of the Eqs. (0) then reduce to the following two conditions: 

ip K (2k 2 + k' 2 )\ , s 

30b*-*) J P = °' (12a) 

^ + 2^) P = ' ( 12b ) 

which give the wave number fe' and the polarization amplitude P. 

On the other hand, Eq. (g) for the electric field refers to a one-atom quantity P + = Pi. There is no evident reason 
to drop the delta function contribution in G, so we use Eq. (Q) as it is. The total electric field is assumed to be of 
the form Eee tk ' z . With the choice (Egg), the vacuum type contributions e ikz indeed cancel, and the polarization 



amplitude is related to the amplitude of electric field by egE = k P/(k — k ). With the help of Eq. (12a) we find 



71 - 1= (1-W3) ' (3) 

where the refractive index n is defined in a familiar manner, k' — nk. Because we are dealing with a linear theory, 
the electric field and the polarization are naturally related by P + = eoxE + ; but in addition we now find the familiar 
relation x = ra 2 — 1 between susceptibility and refractive index. 

It should be noted that keeping or ignoring the delta function of G in Eq. (|^) simply makes a difference between 
two different definitions of electric field inside a dielectric medium. This is a thorny and often an inconsequential 
issue. The final measurements are usually carried out outside the dielectric, whereupon the delta function does not 
contribute anyway. Our choice has the advantage that it aligns with the standard conventions of electrodynamics. 

While we earlier [Q obtained the column density result, n 2 — 1 = \ — P a i Eq. ( |l3| ) is the classic column density 
result with the added LL local-field correction |l4|]. As far as our microscopic argument is concerned, the correction 
emerges as a result of divergent level shifts of the atoms due to dipole-dipole interactions. Within the present approach 
the main approximations for the LL shift are two: collective linewidths and line shifts are ignored, and correlation 
functions for ground state atoms as appropriate for a BEC are assumed. Atom-atom correlation functions, though 
trivial, are formally taken into account to infinite order. In comparison, Ref. ||ll[| strives at taking into account 
collective linewidths and line shifts, but in exchange makes approximations concerning atomic correlation functions. 

Apart from heralding the physics insights that originally lead to the Lorentz-Lorenz formula, our results should 
have practical consequences in the ongoing studies of the condensates. The LL shift could, and possibly should, be 
studied as the first correction to the column density arguments, which so far have been the exclusive tool in the 
analysis of the experimental results. As another example, we conjectured in that the entire linear hierarchy of 
Eqs. (0) for I = 1,2,..., N may be solved by solving numerically the classical equations for the coupled system of 
electromagnetic fields and N charged harmonic oscillators. The position correlation functions p\ could be taken into 
account by repeating the solution for a number of initial configurations chosen from the proper stochastic ensemble, 
and averaging the results. This type of a Monte-Carlo approach would seem difficult in the presence of delta functions, 
because they dwell in a set of measure zero in phase space of the positions of the atoms and yet may substantially 
affect the results. However, as the troubling delta functions have now proven inconsequential, the door is ajar for an 
essentially exact numerical solution of the optical response of a dense, near-resonance BEC. 
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